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High performance control design for a flexible space structure is challenging since high 
fidelity plant models are difficult to obtain a priori. Uncertainty in the control design 
models typically require a very robust, low performance control design which must be tuned 
on-orbit to achieve the required performance. Closed loop system identification is often 
required to obtain a multivariable open loop plant model based on closed-loop response 
data. In order to provide an accurate initial plant model to guarantee convergence for 
standard local optimization methods, this paper presents a global parameter optimization 
method using genetic algorithms. A minimal representation of the state space dynamics is 
employed to mitigate the non-uniqueness and over-parameterization of general state space 
realizations. This control-relevant system identification procedure stresses the joint nature 
of the system identification and control design problem by seeking to obtain a model that 
minimizes the difference between the predicted and actual closed-loop performance. 


I. Introduction 

Implicit in the design of high performance control systems is the availability of an accurate model of 
the system to be controlled. System identification methods are typically directed toward obtaining accu- 
rate nominal models of an open loop system with bounds on the associated uncertainty. Recognizing this 
dependence on accurate nominal design models, robust control theory has been developed to accommodate 
modeling errors. Recently attention has been drawn to the fact that the issues of system identification and 
control design must be treated as mutually dependent. 1 Since the magnitude and spectrum of excitation 
forces are controller dependent, an appropriate model for control design cannot be determined independent 
of the controller. Since the model that is most appropriate depends on the control design, the open- loop 
response of a model is not sufficient to indicate the fidelity of the model for control design. A better approach 
from a closed loop performance perspective is closed-loop system identification, that is the identification of 
the open-loop plant given closed-loop response data and knowledge of the compensator dynamics. 

An approach to iterative closed-loop system identification and controller redesign is given by Liu and 
Skelton where a state space model is estimated using the q-Markov Cover algorithm and a controller is 
designed using the Output Variance Constraint algorithm. 2 q-Markov Cover theory describes all realizations 
of a linear system which match the first q Markov parameters and covariance parameters of the true system 
generated by pulse responses. First the entire closed- loop system is estimated and then using knowledge of the 
compensator dynamics, the plant is extracted. However, the identified open-loop plant has dimension equal 
to the controller dimension plus the dimension of the identified closed- loop plant. Model order reduction is 
used to remove the superfluous states. A similar approach is to formulate a method which first obtains the 
closed- loop Markov parameters and then the open- loop Markov parameters are recovered. 3 A discrete time 
state space realization is then obtained from the Markov parameters. 

A classical approach to parameter estimation which has been extended to closed-loop system identification 
is the prediction error approach. 4 This optimization method estimates the parameters of a linear system by 
minimizing the squared sum of the errors between the actual measurements and the predicted measurements. 
Recent work extended the prediction error method to closed-loop identification for multivariable systems. 
Building on a method developed for estimation of the parameters of an open-loop system in canonical form 
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from open-loop data, a new procedure for closed-loop system identification was developed and demonstrated 
therein. As with all local optimization methods that utilize inverted Hessians and gradients for searching, 
the optimization method was found to be severely limited by the accuracy of the initial guess. In this paper, 
a global optimization step utilizing genetic algorithms is introduced to determine a coarse optimal solution 
to initialize the local optimization process previously developed. 5 

II. Closed Loop System Identification 

The traditional approach to control design is to obtain a nominal model of the open loop plant, P, which is 
the basis for control design. Since P is an approximation of the true plant, P, the model based compensator, 
Cp, must provide a certain level of robust stability (i.e., Cp must internally stabilize P and P). In the 
system identification process an attempt is made to bound the model error || P — P ||, which determines the 
amount of robustness required by the control design. It is also desirable that the controller provide some 
level of performance robustness; that is the achievable performance should not differ significantly from the 
nominal design performance. 

An upper bound on achievable performance is 6 

II J{P,Cp) || < 

\\J(P,C P )\\ + \\J(P,Cp)-J(P,Cp)\\ (1) 

where 

|| J(P,Cp) || : the achieved performance 

|| J(P, Cp) || : the nominal performance, and 

|| J(P,Cp) — J(P,Cp) || : the performance differential. 

The choice of a specific performance metric J and norm || * || is determined by the control design methodology. 

High performance control design for a flexible space structure is especially challenging since high fidelity 
nominal plant models are difficult to obtain. To achieve high performance requires: 

• high nominal performance 

(•HI J{P,Cp) || small) 

• robust performance 

(HI AP,Cp)~j(P,Cp) ||«|| j(p, Cp) ||) 

The second condition is actually a requirement on model fidelity and indicates that the nominal closed-loop 
system model must closely approximate the actual closed-loop system performance when the compensator 
Cp is used. Therefore, the fitness of the nominal model is a function of the compensator and must be 
judged from a closed-loop perspective. Interestingly, good fitness is not guaranteed by good open- loop 
model matching nor is it precluded by poor open-loop model matching. 1 

The closed-loop system identification process may be considered as one step of an iterative closed-loop 
system identification and control design procedure. It is assumed that a moderately accurate dynamic model 
of the system to be controlled is available for the initial low authority controller design. However, this initial 
design model is not of sufficient fidelity to permit high authority control design. The objective of the closed- 
loop system identification procedure is to refine the initial control design model based on closed-loop response 
data. 

In the development of a closed-loop system identification method, consideration must be given to the 
non-uniqueness of the state space realization of the system dynamics. Although there are an infinite number 
of equivalent state space realizations for a system, a system with n states, nu inputs, and ny outputs can be 
uniquely expressed with a minimum of n(nu + ny) parameters. In order to circumvent the non-uniqueness 
problem, this paper will implement an approach which realizes the open-loop system matrices in a unique, 
minimal form and directly identifies the canonical parameters from closed-loop response data. 

Proper selection of the objectives of system identification and control design further stresses the joint 
nature of the identification and control problem. The objective of the closed-loop identification procedure 
employed herein is to obtain a model P that minimizes the performance differential || T — T || 2 where T 
is the actual closed-loop system and T is the identified closed-loop system. This system norm cannot be 
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evaluated since T is not known. However the actual and predicted closed-loop measurements are known and 
an equivalent objective is to minimize the prediction error of the closed-loop system, || y — y || 2 - The actual 
and predicted closed-loop system outputs, y and y, respectively, are determined for the same set of inputs. 
The least squares cost functional for control dependent closed-loop system identification then is 

J = \ J (y — y) T W (y — y)dt (2) 

where W is a constant matrix chosen to weight the relative importance of different measurement outputs. 
The control objective is matched with the identification objective by designing the controller to minimize 
the H 2 criterion || T || 2 - 


III. State Space Canonical Realization 


A state space representation of the open loop plant may be given by: 

(3) 

(4) 


z = Tz + Qu , 2(0) = zo 
y = Hz 


The objective is to identify some IF, Q , H, and z 0 which represents the dynamics of the unknown system based 
on knowledge of the inputs, u(t), and noisy measurements, y(t). An estimate of the unknown parameters 
may be obtained by minimizing the cost functional given in Eq. 2. Directly minimizing J for the general 
state space realization results in a nonlinear optimization process, but in the absence of measurement noise, 
a linear formulation may be obtained. Recognizing that for a perfect model the output in Eq. 4 will exactly 
equal the measurements and y — y = 0, this difference can be fed back to the model with arbitrary gains K, 
and M according to 


which can also be written 


by use of the definitions 


i 

= Tz + Gu + JC(y — Hz) 

(5) 

y 

= Hz + M{y — Hz) 

(6) 

z 

= Hnz + Gnu + 5Qu + ICy 

(7) 

Z{ 0) 

= zno + 5z 0 

(8) 

y 

= HNZ + My 

(9) 


T n = T - JCH 

(10) 


Gn = G — SG 

(11) 


H n = (/ — M)H 

(12) 


zno = z 0 - 5z 0 

(13) 


The elements of (T — K.H) and (J — A4)H may be chosen independently of the unknown elements of T and 
H by using a maximum of n * ny parameters in /C and A4 when the structure of the system is in a specific 
form. As a consequence, Tn, Gn , Hn , and zno are chosen and the unknown parameters are contained in 
/C, A4, SQ , and Szq. The vector of uncertain parameters to be estimated is 


7 = 


vec(IC) 

vec(SG) 

vec(M) 

vec(Szo) 


(14) 


A canonical realization of the system dynamics is necessary so that the unknown parameters in T are 
coefficients of measured states. This paper utilizes the observer canonical form. 7 It can be shown that if 
the plant dimension is an even multiple of the number of outputs and the first n rows of the observability 
matrix are linearly independent, then the realization is canonical and H will consist only of ones and zeros. 
Otherwise, some elements of H will be included as unknown parameters in the estimation procedure. 
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A. Closed-Loop System Realization 

A similar approach may be taken for identification of the closed-loop system parameters by expressing the 
plant in observer canonical form and exploiting the structure of the closed- loop matrices. For the plant given 
by 


x = Ax + B\w + B 2 u 
y = C 2 x 


and a dynamic compensator 


x c — A c x c -I - B c y 
u = —C c x c 

the resulting closed-loop dynamics are 


X 


A 

— B 2 C c 


X 

X c 


B c C 2 

A c 


X c 


B\ B 2 


W 

0 0 


u 


1 

0 

1 


X 

1 — 1 

0 

1 


X c 


(15) 

(16) 


(17) 

(18) 


(19) 

(20) 


If the plant (A, B 2 , C 2 ) matrices are expressed in observer form, then the C 2 matrix consists of ones and 
zeros and the unknown parameters in A are coefficients of the measured transformed states. (As stated 
earlier, in some cases the transformation may not be canonical resulting in the C 2 matrix having additional 
nonzero elements, but these parameters can be estimated as well.) Comparison of Eqs. 3-4 with Eqs. 19 - 
20 indicates that 


A ~B 2 C c 
B c C 2 a c 


(21) 


g 


B\ B 2 

0 0 


(22) 


n = 


c 2 0 

0 1 


(23) 


Since the only unknowns in T and Q are the plant matrices A, B\, and B 2 , the unknown parameter 
matrices are defined as 


AC 


ACn -IC 12 C c 
0 0 


(24) 


65 


5 Q 11 6 Q 12 

0 0 


(25) 


In the case of a true canonical form for the plant, 'H is completely known. Compatibility of AC and Ti 
with T requires TL to have the form given by Eq. 23 which implies that the compensator state vector 
must be included in the measurement vector used in the identification process. Note that ACn corresponds 
to the unknown parameters in A, AC 12 corresponds to the unknown parameters in B 2 , SQu corresponds to 
the unknown parameters in B\, and 6Q \ 2 = AC \ 2 . Hence for closed-loop system identification, the unknown 
parameter vector becomes: 


7 = 


uec(ACn) 
uec(AC 12 ) 
vec(6G n) 
vec(6z 0 ) 


( 26 ) 
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Local optimization methods have been shown to be quite effective for parameter estimation when a good 
initial estimate is available . 5 * However in the absence of good initial estimates, these methods tend to be 
unsuitable due to numerical sensitivity. Local optimization methods often require the inversion of a large 
data matrix which tends to be ill-conditioned. Ensuring full-rank of the data matrix requires sufficiently rich 
excitation in all of the modes to be identified which is often quite difficult in practice due to limitations such as 
sensor and actuator dynamics. Methods which do not involve matrix inversion such as genetic algorithms have 
potential for the closed-loop parameter estimation. An additional benefit of genetic algorithms is that the 
system identification can be operating in the background as part of an autonomous identification/controller 
tuning process. Hence this paper presents a genetic algorithm for the initial (coarse, global) estimation step 
with the objective of providing a sufficiently accurate initial guess for the final (local) optimization step using 
the equation error method. 


IV. Genetic Algorithms 

Scientists and engineers often borrow concepts from nature to motivate and inform the design process. 
One such innovation is the computational field of genetic algorithms. Genetic algorithms have seen increased 
attention over the last two decades in a wide range of engineering applications from aircraft design to 
operations research. Genetic algorithms are mathematical methods for global search and optimization that 
are based on the mechanics of natural selection and the science of genetics . 8 

In nature, natural selection operates on biological systems to determine the fittest species for a given 
ecological niche. The fittest species tend to survive and propagate the genetic information that rendered the 
species most fit. In this sense, natural selection has a stabilizing effect that tends to increase the presence of 
the fittest species through subsequent generations by selecting the fittest members for propagation. Mating 
within the population also tends to mix the genetic information of the population such that beneficial traits 
of parents can be combined to improve the fitness of the offspring. Hence the operations of selection and 
reproduction tend to preserve the fittest members of a species and in some instances improve the fitness 
of individuals within the population through genetic mixing. However selection and reproduction both 
operate on a fixed gene pool with fixed genetic information. Truly innovative biological designs require the 
introduction of new biological information. 

The only known natural mechanism for introducing new biological information is mutation, which es- 
sentially modifies the genetic code of a species. Mutations are random, relatively rare, and almost always 
harmful to the species. However the exceptional beneficial mutation introduces new genetic information 
to the gene pool which may lead to more fit members to be selected and reproduced. This paradigm of 
biological optimization has been fruitfully applied to engineering optimization methods in the general field 
of evolutionary computation and genetic algorithms in particular. 

At the cellular level, biological information is coded in one or more chromosomes which serve to describe 
the function and construction of the biological system. Chromosomes are comprised of strings of DNA 
which can be divided into substrings called genes. Each gene is a functional block of DNA that encodes the 
biochemical workhorse of the cell - protein molecules. The genetic code is simply a mapping from DNA 
strings to individual proteins and the genome is an organism’s set of genes that comprises the set of biological 
characteristics or traits of the organism. So a trait (such as hair color) is encoded by a gene and the different 
possible realizations of that trait (such as brown, black, etc.) are called alleles. 

To carry over the biological paradigm as a computational paradigm, the chromosome may be thought of as 
a candidate solution to the optimization problem. The chromosome is a string formed by the concatenation 
of individual genes, or parameters to be optimized. Alleles correspond to the values that may be taken by the 
parameter (either one or zero in a binary representation, a real number in a floating point representation). 

Genetic algorithms consist of coding a population of candidate solutions and then operating on the 
population to determine a new population for the next generation. Several generations are typically iterated 
upon until the population converges to its most fit set of solutions. The operators for one generation typically 
include selection, recombination, and mutation. 

Selection implements the operation of natural selection by evaluating the “fitness” of each member (from 
the cost functional to be optimized) and giving more copies of the better solutions in the population for the 
next generation. Selection takes many forms such as roulette wheel selection, where the probability of an 
individual being propagated to the next generation is weighted by the fitness; tournament selection, where 
the random members compete “head to head;” and truncation, where the top “n” members are selected and 
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copied multiple times to fill the population for the next generation. 

Recombination implements biological reproduction by combining two parents to form two children for the 
next generation. Two strings chosen at random are cut at a random location in a gene (for binary encoded 
strings) or chromosome (for floating point encoded strings) and the remaining portions of the strings are 
swapped. Also known as crossover, this step serves to randomly shuffle the information in the gene pool, 
but it does not introduce new information in the gene pool. For example, if every member in a population 
has a “1” in the i th string location, then neither reproduction or recombination will produce a “0” in that 
location. 

The introduction of new information is accomplished through the third operator, mutation. Mutation in 
a genetic algorithm is accomplished by randomly altering an allele. For a binary encoding, mutation simply 
flips a bit at random locations as a function of the probability of mutation. For a floating point encoded 
string, mutation perturbs a random allele is some fashion such as a random variation over some range of 
parameter variation. Taken together, these three operators serve to advance the population from generation 
to generation in a manner that tends to improve the fitness of either the population as a whole or the best 
members of the population. 

The next section presents an example application of genetic algorithms to the problem of parameter 
estimation for closed loop system identification. 


V. System Identification Implementation 


Standard local optimization methods typically involve differentiation of an objective function and inver- 
sion of matrices that may tend toward ill-conditioning. Genetic algorithms are much simpler to implement 
than standard methods because genetic algorithms merely evaluate the objective function. This distinction 
also implies that genetic algorithms are better suited for global searches since gradients inherently provide 
only local topological information. Clearly the over-parameterization of a general state-space representation 
will introduce difficulties for the genetic algorithm search just as it does for standard optimization methods. 
Implementing the canonical realization for the closed loop system identification mitigates this difficulty and 
permits the genetic algorithm to search for the best set of unknown parameters for the minimal realization. 
Hence the objective of the genetic algorithm for closed loop system identification is to determine the plant 
parameters which minimizes the error between the predicted and measured outputs as quantified by the 
objective function in Eq. 2. 

Given an initial closed loop plant described by To, Qo, and Ho, the updated plant matrix T new may be 
given by 

T n ew = To - AT (27) 


Since the unknown plant parameters to be estimated are elements of (A, Bi, B 2 ), inspection of Eqs. 21 and 
10 indicates that 


AA -A B 2 C c 
0 0 


= KH 


(28) 


With K, and H compatibly partitioned as defined in Eqs. 24 and 23, Eq. 28 yields AA = JCnC 2 , or 


K n = A AC+ 


(29) 


where the right inverse (pseudo-inverse) of C 2 is given by 

c+ = cj {ciczy 1 


(30) 


The variation in A can be treated as a multiplicative percent variation from the initial matrix Ao by the 
definition 


AA — A^ Aq 


(31) 


where A^ is the matrix 


A^i = diag [8Ai\ , Vi = 1 : n * p 


(32) 


and the ” diag” notation indicates a diagonal matrix with the elements in brackets along the diagonal. Each 
8 A-i represents the percent uncertainty in the i th uncertain element of the open loop A matrix. Using the 
canonical realization in this manner allows us to isolate the individual uncertain elements of the open loop 
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plant matrices from within the closed loop realizations as linear multiplicative coefficients. Likewise, Eq. 28 
yields 

A B 2 C c = IC 12 C c H 22 (33) 

For this case where the realization is canonical, 'H 22 is the identity matrix and hence 

/C 12 = AB 2 = A b 2 B 2 q (34) 

where B 2 q is the initial value of B 2 and A b 2 is a diagonal matrix of percent variations analogous to Eq. 32. 
Finally, the updated Q plant matrix is given by Eq. 11 which by direct inspection yields 

SG U = ABi = A Bl B h o (35) 

where B\ q is the initial value of B\ and A Bl is a diagonal matrix of percent variations analogous to Eq. 32. 

This representation allows the chromosome to be a concatenation of the percent variations of the indi- 
vidual uncertain plant parameters. Beginning with a population of chromosomes, each binary chromosome 
string is decomposed into substrings, or genes, that represent the percent variation of each plant parameter. 
The substrings are decoded into floating point numbers that are the diagonal elements of Aa,A B i , and 
A s 2 . These floating point numbers are then used in Eqs. 29 - 35 to form the parameter vector 7 defined 
by Eq. 26 (neglecting uncertainty in the initial conditions Sz 0 ). Hence each chromosome in the population 
is mapped to a corresponding 7 from which the closed loop plant matrices are formed. A time response of 
the corresponding closed loop plant is compared to the measured time response of the true system and the 
fitness of each chromosome in the population is finally determined by Eq. 2. 

VI. Application Example 

The four disk structural system 9 is an excellent example of a flexible space structure and is used herein 
as a test case for closed loop system identification using genetic algorithms. The four disk model can be 
configured with a rigid body mode and coupled nonminimum phase flexible modes (as a consequence of 
noncolocated sensors and actuators). As illustrated in Figure 1, four disks are rigidly attached to a flexible 


Ml, 01 M2, 02 M3, 03 M4, 04 



12 3 4 


Figure 1. Four Disk System 


axial shaft with control torque applied to selected disks and the angular displacement of selected disks 
measured. The equations of motion may be written as 


or 


hOi + Ki(0\ — 0 2 ) = Mi 

I 2 § 2 + K 2 {9 2 - 6*3 ) - K x (6>i - d 2 ) = M 2 

+ K 2 (93 — 6 b) — K 2 (9 2 — 63 ) = M 3 

iJi- K 3 {9 3 - 9 4 ) = M 4 


Mq + Kq = Bu 


(36) 


(37) 
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where the mass matrix is 


h o 0 0 ■ 

0 I 2 0 0 

0 0 h o 

0 0 0 h _ 

and the stiffness matrix with ( GJ/L)i = A) is given by 

Ad -A' 2 0 0 

K -Ad Ad + Ad -Ad 0 

0 -AT 2 Ad + K 3 -K 3 

.0 0 -K 3 Ad 



(38) 


(39) 


The generalized displacements are the angular displacements of the disks, q T = [0\ 0 2 0 3 64 ], and the input 
vector consists of the moments applied to each disk, u T = [M\ M 2 M 3 M 4 ]. Defining the state vector as 
x T = [q T q T ] results in the state space formulation 


x = Ax + B 2 u 


(40) 


where 


A = 


0 / 
M~ 1 K 0 


B 2 


0 

M~ l B 


(41) 


Note that Eq. 19 includes an explicit input coefficient matrix, B\, for the disturbance input, w. This 
disturbance input typically consists of noise or uncertainties at the plant input and is used as a design 
parameter for controller synthesis. In the control design context, the B\ matrix includes the control design 
weights or states associated with noise or uncertainty models and does not explicitly model the dynamics of 
the flight system. Moreover, rarely is a spacecraft instrumented to directly measure the disturbances even 
if the disturbance sources are fully known. Since the disturbance time histories are not measured (required 
for system identification), B\ is not included above in Eq. 40. If on-orbit dynamics characterization is 
anticipated, then excitation actuators (disturbances) may be included in the system design. But these 
actuators can also be used for feedback control as well and as such would be accounted for in the B 2 matrix. 
Hence, all measured inputs are accounted for in the B 2 matrix and the estimation of B 1 is not required. 
This has the beneficial effect of reducing the number of parameters to estimate and hence reducing the 
chromosome length. 

For this example, the true plant is described by an identity mass matrix and the stiffness values are 
K\ = 25, K 2 = 12.5, and K 3 = 1. These values result in an open loop system with undamped modes at 
0.145 Hz, 0.347 Hz, 0.688 Hz, and 1.156 Hz. A noncollocated sensor/actuator configuration is adopted with 
the outputs consisting of the angular displacement of disks three and four while the control torque is applied 
to disk one. 

Case A 

The first case begins with an estimated model that has a 10% error in the stiffness matrix and a 5% 
error in the mass matrix. With eight states, one input, and two outputs, the chromosome was a string of 
n* (p + m) = 24 parameters encoded as an eight bit integeter. The eight bits mapped into a range of ±10% 
for each parameter. 

An initial population of 4000 chromosomes was generated with random binary bit strings. In order to 
preserve the best member of each generation, an elitist strategy was implemented where the most fit 10% 
of the population was directly copied into the subsequent generation. Tournament selection was used with 
random pairing of two chromosomes; the probability of crossover was 0.75, and binary mutation occurred 
for 50% of the parameters (the probability per bit was 0.0625 with eight bits per parameter). This was an 
atypically high probability of mutation because the objective was not to converge the entire population to 
the most fit member as in the typical genetic algorithm. For this application, the one most fit member was 
sought and a high mutation rate was employed to continually infuse new information content into the gene 
pool. 

Feasibility of solutions was enforced by checking the closed loop stability for each member of the popula- 
tion and rejecting those members that rendered the closed loop unstable. Clearly the difficulty of populating 
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the gene pool with feasible plants is a function of the control system. For this reason also a robust controller is 
most suitable for closed loop system identification (as described earlier). A low authority LQG controller was 
designed for the initial plant and held fixed for the plant estimation phase. In a flight application, the best 
plant identified during this step would be used for a controller update. This iterative identification-tuning 
process could be repeated until acceptable performance was achieved. 

Figure 2 presents the evolution of the gene pool 
for 350 generations for two runs. Each run begins 
with the same initial erroneous closed loop plant but 
with different random initial populations and differ- 
ent disturbance time histories excited through the 
control actuator. The forty second disturbance in- 
put time history is a sum of four random time histo- 
ries (with intensities of 0.25, 0.5, 0.75, and 2). The 
disturbance input is band limited by a 6th order but- 
terworth filter with a break frequency at one-half the 
sample frequency of 10 Hz. 

Several different criteria can be applied to deter- 
mine the suitability of an estimated plant. Fitness 
alone is difficult to correlate to absolute model er- 
rors because the cost is an integral of the output 
error squared over the duration of the time history. 

Nor is fitness is an exact basis for comparison be- 
tween models since different time histories are used 
for different runs. However, because the dynamics 
are realized in a canonical set of basis vectors, the 
parameter values themselves are equivalent for equivalent models. But determining the significance of an 
error in those parameters is not straightforward. Table 1 presents the parameters the true plant, the initial 
plant, and the best estimate for two runs (370 generations and 800 generations, respectively). The first eight 
rows are elements of the A matrix and the final three are elements of the B i matrix. Note that in the first 
case, the best estimate in terms of fitness was achieved with larger errors in two elements of the A matrix. 
In the second case, errors in all elements of the A matrix were improved upon. This is an indirect indication 
that not every element of the A matrix contributes equally to the cost and larger error in some parameters 
is more acceptable than others. Note that errors in the elements of B\ tend to have less effect on the cost 
than errors in A. 


True 

Init. 

Est. 1 

Est. 2 

Init. % error 

Est. 1 % error 

Est. 2 % error 

-625.00 

-685.94 

-639.40 

-629.72 

9.75 

2.30 

0.76 

-62.50 

-65.48 

-62.32 

-62.37 

4.76 

-0.29 

-0.21 

25.00 

27.44 

25.55 

24.80 

9.75 

2.22 

-0.79 

1.00 

1.05 

1.09 

1.02 

4.76 

8.58 

2.26 

156.25 

171.49 

169.27 

163.62 

9.75 

8.33 

4.71 

156.25 

163.69 

168.25 

162.73 

4.76 

7.68 

4.15 

-12.50 

-13.72 

-12.98 

-12.66 

9.75 

3.85 

1.27 

-14.50 

-15.19 

-14.86 

-14.59 

4.76 

2.50 

0.61 

37.50 

37.42 

37.66 

37.22 

-0.23 

0.44 

-0.74 

1.00 

0.95 

1.01 

1.01 

-4.76 

1.03 

-0.39 

-1.00 

-1.00 

-1.01 

-0.97 

-0.23 

0.67 

-3.24 



Figure 2. Fitness Evolution for Case A 


Table 1. True, Initial, and Estimated Parameters for Case A 


Another means of comparison is the eigenvalues of the open loop plants. The eigenvalues have a more 
direct effect on the time response of the system than do the individual elements of the plant matrices and as 
such are a valuable basis for comparison. Figure 3 shows the very significant reduction in error associated 
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with frequencies accomplished by the genetic algorithm parameter optimization. This demonstrates the not 
surprising fact that like the individual matrix elements, not every mode contributes equally to the cost. As 
with the parameter errors, second run yielded a more accurate model by this means of evaluation. Both 
estimates are quite accurate in terms of parameter values and open loop eigenvalues for this case. 

Case B 

The second case begins with an estimated model that has a 25% error in the 
stiffness matrix and employs a ten bits per parameter binary encoding mapped 
into a range of ±50% for each parameter. Otherwise the parameters are the 
same as Case A. Table 2 presents the comparison of parameter estimates which 
indicates a significant reduction in error with most parameters, although some 
parameter errors are actually increased. Each of the four modes in the initial 
model erred by 9.0% while the estimate modes erred by -1.0436%, 0.50%, .027%, 
and -0.38%. This case clearly demonstrates how the mode frequencies of the 
model can very accurately match the true plant even though the errors in some 
evidently insignificant parameters can be quite large. 


True 

Init. 

Est. 

Init. % error 

Est. % error 

-625.00 

-885.77 

-545.06 

41.72 

-12.79 

-62.50 

-74.40 

-58.29 

19.05 

-6.73 

25.00 

35.43 

35.07 

41.72 

40.27 

1.00 

1.19 

1.78 

19.05 

78.46 

156.25 

221.44 

136.48 

41.72 

-12.65 

156.25 

186.01 

140.10 

19.05 

-10.34 

-12.50 

-17.72 

-15.82 

41.72 

26.55 

-14.50 

-17.26 

-18.35 

19.05 

26.55 

37.50 

42.52 

32.98 

13.38 

-12.06 

1.00 

0.95 

1.00 

-4.77 

-0.43 

-1.00 

-1.13 

-1.57 

13.38 

56.55 



Figure 3. Error in Open 
Loop Plant Eigenfrequen- 
cies for Case A 


Table 2. True, Initial, and Estimated Parameters for Case B 


Case C 

The final case includes a nonzero and uncertain 
damping. This case begins with an estimated model 
that has a 35% error in the stiffness matrix and 50% 
error in the damping ratio with the binary strings 
mapped into a range of ±75% for each parameter. 
Otherwise, the parameters are the same as Case B. 
Table 3 presents the comparison of parameter esti- 
mates while Figure 4 presents a comparison of eigen- 
values for Case B. Note that only the A matrix pa- 
rameters for the first two genetic algorithm runs are 
included while two additional runs are included in 
the mode frequency error bar plot. The first esti- 
mate is the most accurate of the four runs consid- 
ered, each of which demonstrated considerable er- 
ror in the frequency of the fourth mode. This case 
shows the difficulty encountered with significantly 
larger uncertainties in the model. However, stan- 
dard local optimization routines are typically much 
less robust to modeling errors. 5 Future research will 
compare the local optimization alone to the genetic 
algorithm/local hybrid optimization for this example. 



Estimate 4 


Figure 4. Error in Open Loop Plant Eigenfrequencies 
for Case C 
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True 

Init. 

Est. 1 

Est. 2 

Init. % error 

Est. 1 % error 

Est. 2 % error 

-625 

-1033.2 

-882.43 

-947.57 

65.306 

41.189 

51.611 

-68.75 

-93.112 

-138.64 

-148.2 

35.436 

101.66 

115.57 

25 

41.327 

29.48 

34.934 

65.306 

17.92 

39.735 

1.25 

1.7959 

1.0125 

1.0178 

43.673 

-19.00 

-18.578 

156.25 

258.29 

186.14 

233.11 

65.306 

19.132 

49.188 

157.81 

204.08 

147.97 

191.96 

29.319 

-6.2342 

21.64 

-12.5 

-20.663 

-16.497 

-21.678 

65.306 

31.978 

73.426 

-14.625 

-18.898 

-17.609 

-22.569 

29.217 

20.407 

54.321 


Table 3. True, Initial, and Estimated A Matrix Parameters for Case C 
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